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The variance of diflference of photocounts is an established measure of quantum correlations for quantum states 
of light. It enables us to discriminate between the classical correlation of a two-mode coherent state and the 
quantum correlation of a twin-beam state. We study the effect of loss and saturation of the photon-number- 
resolving detector on the measurement of the variance of difference of photocounts. An analytic function is 
derived for this variance, both for the coherent and the twin-beam states. It is found that the variance of 
difference of photocounts is no longer a reliable entanglement measure in the non-linear regime of the detector 
response but it remains useful in some range of values of average photon numbers of the incident light. We also 
quantify the linear regime of the detector with saturation which will be useful for calibration of the detector 
quantum efficiency. 
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1. Introduction 

Generation and characterization of multiphoton entangled states represent one of the prerog- 
atives of modern quantum optics experiments. Highly-entangled multiphoton states such as 
NOON states [1] and photon-number-correlated twin-beam states [2] are promising resources for 
practical applications in quantum lithography [3J, quantum metrology, [1] and quantum cryp- 
tography [5j. Measurement of the degree of quantum correlations accurately is an essential task 
for characterizing states of these types. The recent emergence of photon- number-resolving detec- 
tors (PNRDs) |6l[7j, such as visible-light photon counters (VLPCs) [8], transition-edge sensors 
(TESs) [9j, time- multiplexed detectors (TMDs) [TO] and multipixel photon counters (MPPCs) 
[TTj, makes the direct characterization of multiphoton states by photocounting possible. 

A measurement using a A^-photon-resolving detector can be modeled by a positive-operator 
valued measure (POVM) with + 1 outcomes: {Ho, Hi, . . . , Hat}, satisfying completeness, 
Si^o ~ ^^'^ positivity, Ilj > 0, V « = 0, . . . , iV. The outcomes of the POVM are diag- 
onal in the Fock-state basis and hence phase insensitive. The probability of detecting i photons 
is p'^ = tr(pnj), where the statistical operator p describes the photon state, and knowledge of 
the p[s enables us to reconstruct the photon statistics. For a perfect detector with no loss (unit 
quantum efficiency) and no saturation, {Hi = \i) = 1,2,...}. Then the probabilities p[ 
should reflect the true frequency of occurrence of i photons received, pi, in the limit of a large 
number of detection trials. However, because of loss and noise in the detector the measurement 
is a proper POVM rather than a von Neumann measurement and the measured probability p'- 
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is a mixture of the actual photon statistics, 

Pi = Y^ WijPj , (1) 

3 

where Wij is the conditional probability of detecting i photons when j photons impinged on 
the detector. The conditional probabilities model the measurement process in the detector; in 
the POVM picture, they are related to the outcomes by Wij = (j| Ilj where \j) is the ket 
of the Fock state with j photons. The set of Wij has been used to quantify the capability of 
photon number discrimination of the PNRD |12j . and the objective of detector tomography is 
to determine the conditional probabilities Wij in a model-free way |13H15j . However, detector 
tomography relies on the ability to prepare a tomographically complete set of states and com- 
plicated optimization techniques. Therefore, parameter estimation for instrument calibration 
requires suitable models for the POVM |1 H [T6l - [20] . 

Highly entangled states can be produced in various nonlinear optical processes [21]. At high 
intensities, the photodetector response becomes nonlinear due to saturation. For example, in a 
silicon multi-pixel photodetector consisting of a spatial array of avalanche photodiodes (APD), 
saturation of individual APDs in the silicon multi-pixel detector results in the saturation of 
the array. TES is another alternative for PNRD, which exploits bolometric methods for photon 
detection. TESs operate at cryogenic temperatures (below lOOmK) and have high quantum 
efficiency (~ 95%). Their saturation at high intensities is due to the heating of the detectors 
beyond the superconducting state. Hence, there is a clear need to understand the effects of 
saturation so that one can use these PNRDs in experiments with intense non-classical states of 
light. 

In this paper, we assess the feasibility of using lossy PNRDs for measuring quantum correla- 
tions under joint photodetection in the presence of saturation. The figure of merit that we use 
here is the variance of difference of photocounts (VDP), which has been shown to be a good 
measure of entanglement [22] . The dependence of the VDP on quantum efficiencies of photode- 
tectors has been exploited in the method of absolute calibration of quantum efficiency using 
twin-beam light |23H25j . The effects of saturation on this calibration protocol will be discussed 
here too. 

This paper is organized as follows: In Section 2, the joint detection setup is described. The 
POVM modeling of losses and saturation are presented. The first and second moments of the 
photocounts of these POVM for a state with Poissonian statistics are derived. In Section 3, the 
VDP and its moment operators are described. The VDP of two pairs of two-mode states, one with 
classical and another with quantum correlations, are calculated. The difference of these VDP, 
Q, is introduced as a measure of distinguishability between quantum and classical correlations. 
A tractable analytic form is derived for Q. In Section 4, we introduce the noise reduction factor 
(NRF), a normalized form of the VDP, and show how an absolute calibration of two detectors 
can be performed using the NRF in the linear regime of the detector response. The analytic form 
of the photocount calculated in Section 2 gives a quantitative estimate for the range of mean 
photon numbers in which the detector remains linear. In Section 5, the results are summarized 
and the use and limitations of the measurement model are discussed. 



2. Modeling photon detection with loss and saturation 

The joint measurement scheme is as follows: Two modes of radiation with the mode creation 
operators, ai and 02, are incident on two PNRDs. The resulting photocounts of the respective 
PNRDs are analyzed. The joint detection of quantum correlations with lossy photodetectors has 
been looked into in f22]. Here we extend the theory to include the effects of saturation. For the 
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purpose of the analysis, we use the twin-beam (TWB) state defined as 

oo 

\X) = ^bn\n),\n)^ , (2) 

n=0 

where \n)j ,j = 1, 2 is the Fock state with n photons in mode j and bn is the probabihty amplitude 
of the state \X). It follows from the structure of Eqn. ^ that the two modes of the TWB state 
are perfectly correlated in photon numbers. The TWB state can be generated by spontaneous 
parametric downconversion (SPDC) in a nonlinear crystal pumped by a pulsed laser. The two 
beams created by SPDC comprise inherently of multiple frequency modes. In the limit of a large 
number of these modes, \bn\'^ = exp(— n)n"/?7-!, where n is the mean number of photons in either 
the signal or the idler mode \26\ I27j. As a standard for comparison with classical states, we shall 
use the two-mode coherent (TMC) state \a)i {0)2 where n = |ap. Without loss of generality, we 
can assume that there is no phase difference between the two modes in the TMC state because 
photocounting is a phase- insensitive measurement. 
The mth outcome of the POVM of a lossy detector with quantum efficiency rj is given by |28j . 

00 

(??)|n)(n| , (3) 

n=m 

with 

^,n,n(r?) = ry-(l-ryr-™(^^) , (4) 

for m = 0, 1, 2 . . . and Wm,n{'n) '^^s introduced in Eqn. ([T]). For a PNRD, the measured number 
of photons, m, which is an eigenvalue of the photocount, m, has an upper bound at N. When 
m < N, then its POVM is the same as in Eqn. ([s]), but because the PNRD is unable to resolve 
between N and -|- 1 or more photons, we have 

N-l 

= / - J] , (5) 

m=0 

for the Ath outcome. A is the maximum resolvable photocount. When A or more photons 
are incident on the detector, at most A photons will be registered. When A = 1, this POVM 
reduces to that of the ON/OFF photodetector: {IIo, Hi}, where Hq = X]5^o(-'- ~ I"") (""I 
Ui = I-Uo. 

Given the outcomes of a joint measurement of a two-mode state p, the joint probability dis- 
tribution is given by p{mi,m2) = tr{pllmi Hma)- Moments of this distribution, {m\m2) can 
be calculated via the p- moment operator of the photocount, 

N-l 

mP = TvPUm 
m=0 

N-l 00 

= - E E - mnB.^,n{v)\n){n\ . (6) 

m=0 n=m 

We note that since they are the operatorial moments of a POVM, mP ^ fhP [22] unless the 
outcomes 11^ are pairwise orthogonal which is not the case here. The first two moment operators 
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are 

oo 

m = rjn — Cn{N,r]) \n) {n\ , (7) 

n=N+l 

and 

oo oo 

= ry2 + r/(l - - (2iV + 1) ^ C„(iV, ?7)|n)(n| + 2 ^ D„(iV, ?7)|n)(n| , (8) 

n=N+l n=N+2 

where 



n 

+ 2 

with X = r]/{l — 7]) and -F(a, 6, c\x) is Gauss's hypergeometric function. 

Without saturation, the first and second moment operators are mo = r]n and ttiq = rj'^n? + 

7y(l — r7)n respectively. The variance of the photocounts is then a'^{m) = {rn?) — {m)'^ = rf'a'^{n) + 
rj{\ — T]) (n). It is clear from Eqs. ((7])-([9]) that the saturation adds extra summation terms into 
the moment operators. For a state with Poissonian statistics we have 

(m)p = N - [NeN-iivn) - r?neAr_2(r/n)]e"''" , 
where e„(x) is the exponential sum function 

n p. 

^n(-) = E|r- (11) 

k=0 

The dependence of the average photocount of a single detector on the average photon number 
of an impinging Poissonian light is shown in Figure [TJ When n — )• cxd, the average number of 
photocounts tends to the saturation value with the following behavior: 



ny \{N-l)ln \n-^ 

By contrast when n <^ N, the saturation effect is negligible so the expected detected photons 
behave like that without saturation, 

(m)p = r,n + {n') - {j^^, " 0{n')^ . (13) 

When (rjn)^ ^ (A^ + 1)!, then {m)^ = rjn which identifies the linear regime of the detector. The 
POVM in Eqs. describes the average photocount of a detector expected in the asymptotic 

limit of small and large average impinging photon number. 



Measuring quantum correlations using lossy photon-number-resolving detectors with saturation 



5 




Figure 1. A plot of the detected photocount (ffi) with the average number of photons n in the incident beam for the lossy 
detector with various maximum photocount of the detector, A'^ (solid), and quantum efficiency, r] = 0.5. Quantum efficiency 
which is defined to be the ratio of the detected photocount to the incident photon number, i.e. the slope of the detected 
photon curve, tends to at large n. The response of a detector with just loss is also shown (dotted) for comparison. 



3. Variance of Difference 

Given the two photocounts from the joint detection, their difference is 

D = rhi — rh2 , (14) 

with integer eigenvalues d = 0, ±1, ±2, . . .. The POVM of the measurement of the difference 
photocount, Qd, is 

N-\d\ ( ^q+d ® , d>0 

Qd= { ® n, , d = o , (15) 
'^=0 [ Ug n,_d , d<o 

where 11^ is the mth POVM outcome of Eqn. ([3]). The moments of the difference photocount 
distribution are then 

D = dQd = "T-i ~ ^2 ; 

d 

Z)2 = ^ d^Qd = ml + ml- 2mim2 ■ (16) 

d 

The VDP for a given quantum state is a'^{d) = {D^) - {of. Let us find the VDP for the TMC 
and the TWB states. For the TMC state, we get 



al{d) = al{mi) + cr^(m2) , 



(17) 
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where a'^{mj) = {m^) — {mj)\ For the TWB state, the quantum correlations gives nonzero 
terms in the covariance of the joint photostatistics and this shows up in the VDP, 

'^lid) = crlimi) + al{m2) - 2 (mima)^ + 2 (mi)p (m2>p , (18) 

where 



{■mim2) X = Vimri'il + ri) - ^ |6nP??inCn(iV2, ??2) - ^ |6n|^??2nCn(A^i, ?7i; 

n=N2 n=Ni 



oo 



+ \bn\^CniNi,m)CniN2,m) ■ (19) 

n=max{Ni,N2) 

In Figure [2| we show the effects of loss and saturation on measuring the VDP for the TMC and 
the TWB states with different values oi Ni = N2 = N and ?yi = ??2 = V- Two observations can 
be made. First, the saturation causes the variance to decrease with increasing average photon 
number for both the TMC and the TWB states and asymptotically approach zero. This means 
that the VDP cannot provide a reliable discrimination of classical and quantum correlations at 
values of n that are too large. Second, loss degrades the measurement of the quantum correlation 
since with increasing loss, the VDP for the TWB becomes closer to that for the TMC. 

However, there exists a range of n for such a discrimination to be possible. Let us define a 
quantity Q by means of 

Q = Q{Ni,N2,m,V2,n) = al{d)-al{d) . (20) 

Q is an indication of how good the discrimination between the classical and quantum correlations 
is. Figure |3] shows the behavior of Q versus n for the TMC and TWB states. The optimal 
discrimination happens for the largest Q values for some nmax that we can solve for if we know 
the quantum efficiencies, rj, and the maximum photocount of the detector, A^, of the detectors. 
The optimal discrimination does not always occur when rj = 1 due to saturation. When the mean 
photon number is high and saturation is likely to occur, having some loss in the detectors will 
help to set off the effect of saturation. However, loss will degrade the quality of the measurement 
of the VDP. As such, there is an interplay between loss and saturation. 
The analytic forms for the VDP for the TMC and TWB states measured by a PNRD with 



saturation and loss in Eqs. (17)-(19) are the main results of the paper. Using these results, 
we analyzed the VDP for different values of r] and N thus illustrating the significance of loss 
and saturation to the measurement of the VDP. Thus, such an analytic analysis is helpful for 
interpreting experimental data. We will also show in the next section, how we can use the features 
of saturation in the analytic model in the absolute calibration of the PNRDs. 



4. Absolute Calibration of PRNDs using the NRF 

Sometimes it is helpful to consider a related measure of quantum correlation — the noise reduc- 
tion factor (NRF), 
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For the TMC state, NRF = 1 with any value of quantum efficiencies for the detectors. For the 
TWB, NRF = with detectors of perfect quantum efficiencies and NRF = 1 — -^^z^ with lossy 
detectors, where rji and r]2 are their respective quantum efficiencies. This maLes the NRF a 
useful measure of quantum correlations for two- mode states, as well as a calibration measure 
for the quantum efficiencies of detectors [25j . First, we will describe the current experimental 
procedure for absolute calibration which is possible in the linear regime of a PNRD. Then using 
the analytic form for the photocount, a range of values of average photon number for the linear 



regime can be derived. Second, we will suggest the use of the analytic form in Eqn. (17)-(19) for 
absolute calibration beyond the linear regime of the PNRD. 

For the first method, let us consider the scenario in which we are given two detectors with 
unknown quantum efficiencies and saturation characteristics. Absolute calibration, for example, 
can be achieved if we have a TWB state by measuring NRF in the linear regime of the detectors. 
The calibration routine is as follows: 

(1) Measure the photocounts, (mi) and {m2), of the two detectors with increasing n. This can 
be done by increasing the pump power of the laser. From this data, the mean photocounts, 
VDP, and consequently the NRF, can be calculated. 

(2) Let us define the ratio of the two quantum efficiencies, k = In the regime where the 
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mean photon number, n 

Figure 3. (Color online) A contour plot of Q versus mean photon number, n, and loss, r], for a balanced joint detection 
where rji = ri2 = ri and A'l = N2 = 3. In this figure, the maximum Q value is clearly visible. The solid line shows the 
optimal n for a given 77 and the dashed line shows the optimal rj for a given n. The higher the value of Q, the better the 
discrimination between quantum and classical correlations. Note that for a fixed value of n, the best discrimination does 
not always occur at = 1 (zero loss). A higher loss is favored at higher mean photon numbers where saturation gets more 
prevalent to obtain a larger Q value. However, higher losses in the photocount will degrade the measurement of the VDP. 



, IV, UilCll, 

by calculating the ratios of the two photocounts. 



photocounts are approximately linear with respect to n, then, ~ /c. So A; can be found 



(3) For this range of average photon number, NRF = 1 — jri:- Using the measured NRF and 



the value of k found in Step (2) above, r/i can be calculated. Last, we can calculate r]2 



from the values of k and rji found. 

Thus, absolute calibration of the detector is possible in the linear regime of the detector. This 
is quantified by the range of n in which the term linear in n is much larger than the next order 



term in Eqn. (10) 



«1. (22) 



{N + l)\ 

The higher the saturation point is, the broader the linear regime of the detector. Although it 



is difficult to pinpoint exactly how small the value on the left-hand side of Eqn. ( 22 ) must be 



current experiments suggest that values of ~ 0.1 are sufficient [6]. The worst range of values of 
n for the linear regime is given for = 1, where we would require that n < For an APD 
where A" = 1, the quantum efficiency is typically of order 33% \19\ . thus the linear regime is 
where n < 0.6. However, for a VLPC with a typical saturation value of A^ = 10 and quantum 
efficiency rj = 85% [19], we have n < 5.4. So the linear regime of the VLPC is about 9 times 
that of the APD. This puts a quantitative estimate for the linear regime of a PNRD from its 
characteristics which are typical for that type of detector and can be taken from its manual. 
Once an approximate linear regime is determined, a careful calibration can be done in this range 
of average photon numbers. 



In principle, using the analytic forms in Eqs. (10), (17) and (19), we can achieve absolute 
calibration beyond the linear regime of the PNRD with just a TMC state using the following 
protocol: 

(1) Split a pulsed laser beam with a 50-50 non-polarizing beamsplitter into two parts with 
equal intensities to form a TMC state. One has to make sure that the laser produces a 
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coherent state which is free of any additional classical noise, see 

(2) Measure the photocounts of the beam in each arm with a detector where the detection 
intervals are triggered by the laser. Then increase n by increasing the pump power until 
there is no visible change in the mean photocount with increasing n. In this limit, the 
detectors are saturated and we have A'"i = [{mi)^^^] and N2 = [(m2)gg^^], where the subscript 
"sat" indicates values at saturation and [•] denotes the closest integer. 

(3) Calculate the difference of photocounts for each detection interval and then the variance 
of this difference for each n value. 

(4) Using the found values of A^i and we can generate an equation for the NRF of the 
TMC state as measured by the PNRD from Eqs. ([Tt]) and ([l9]). A nonlinear fit of 



the experimental data of the VDP for the TMC state can then be done with the analytic 
form for the NRF using a suitable numerical algorithm. This will yield fitted values of r]i 
and 772 thus completing the calibration of the PNRDs. 

This absolute calibration procedure with the TMC state is possible owing to the features of 
the VDP in the presence of saturation. It would not be possible if detectors operating in the 
linear regime or detectors with just losses are used, because in those cases, the VDP would 
be insensitive to the quantum efficiencies of the detectors. The advantage of using this second 
method instead of the first is that we do not have to rely on a TWB state which are harder 
to generate than a TMC state. The disadvantage is that we have to use a numerical nonlinear 
fit method for the calibration which might yield more than one reasonable fit and the errors 
in these fittings are generally larger than those for a linear regression fit required for the first 
method. 



5. Conclusion 

A model of joint photodetection with lossy detectors considering saturation is presented. Using 
this model, we have derived an analytic form for the VDP which is tractable for calculations. 
We then used this analytic form to study the effects of saturation on measurements of the VDP. 
It has been found that saturation diminishes this measure and the variance of photocounts of 
the TMC and TWB states cannot be distinguished at large mean photon numbers. Loss in the 
PNRDs degrades the quality of the measurement of the quantum correlation. On the other hand, 
loss can offset the effects of saturation at high average photon numbers. Nonetheless, the VDP 
is still a good measure of quantum correlations as long as we are not far from the optimal Q. 

Saturation also limits the range of n in which data is useful for calibration. This has been 
known for some time experimentally, but here we have derived quantitatively what this range is 
in terms of relative values of quantum efficiency and saturation value. With the analytic forms 
for the measured photocounts and VDP, we will also be able to do an absolute calibration with 
a TMC state. The generation of a TMC state is far simpler than that of a TWB state where the 
latter is required for absolute calibration. 

As the focus of this paper is on loss and saturation effects in PNRDs, we have not looked 
at possible extensions of the measurement model. For example, apart from loss and saturation, 
PNRDs can also suffer from dark counts and crosstalk effects, both of which add spurious counts 
to the data. In using the measurement model for accounting for noise in actual experiments, one 
would have to take these effects into account. This can be done by modifying Eqn. ([s]) with the 
appropriate conditional probability matrices representing different noise types. Another useful 
extension of the measurement model is that of a variable saturation value, A^. Here, we have 
considered a fixed saturation value but in practice, a PNRD can have different saturation values 
at different intensities or wavelengths of incident light. A more robust measurement model can 
be obtained by taking this variability into account. 
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